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Abstract 

Computationally efficient algorithms are presented for determining 
single-axis attitude from the measurement of arc lengths and dihedral 
angles. The dependence of these algorithms on the solution of trigono- 
metric equations has been much reduced. Both single-time and batch 
estimators are presented along with the covariance analysis of each 
algorithm, 
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I. Introduction 


Since nearly every spacecraft is spinning during part of its 
life — in particular, at the time of orbit injection — spin-axis attitude* 
estimation is an important segment of almost every mission support 
operation. Indeed, for spin-stabilized spacecraft there is often no need 
(or desire) to determine the complete three-axis attitude at every 
point and, in fact, when accuracy requirements for the spin-axis attitude 
dictate that many measurements taken at different times be processed 
simultaneously, the computation of a three-axis attitude may not even be 
possible. 

Very often, three-axis attitude information is definitive data 
required chiefly by mission scientists and generally processed anytime 
from several days to several months after the receipt of telemetry. The 
need for efficient three-axis attitude estimation algorithms in those 
cases is determined by the definitive data rate. When three-axis 
attitude information is required in real-time for the purpose of attitude 
control, this is usually provided on-board by three-axis gyros (e.g. SMM) 
or on the ground by the spin axis and a third angle, which can be 
obtained by monitoring some other sensor reading such as IR scanner pitch 
(e.g. AEM, Magsat). 

Spin-axis attitudes by contrast are usually required not only as 
definitive data but also by the ground support system in near real-time 
for the purpose of monitoring spacecraft performance and determining 
large scale attitude maneuvers. Thus, the efficiency of a spin-axis 
attitude estimation algorithm becomes a factor in the safety and daily 
operation of the spacecraft. 


* 

Since the single-axis attitude of interest is invariably the spin-axis 
attitude these terms will be used almost interchangeably throughout this 
work. 


16-2 



While a number of highly-efficient algorithms exist for three-axis 
attitude estimation,^ the computation of spin-axis attitude^ is by 
comparison very clumsy. This is largely because the computation of 
three-axis attitude uses complete vector measurements in general and can 
take advantage of the linear properties of Euclidean three-space. The 
computation of spin-axis attitude, on the other hand, must rely on 
incomplete vector information (the measurement of arc lengths and 
dihedral angles) to determine a quantity (the spin-axis) which is 
restricted to the surface of a sphere. Thus, while three-axis attitude 
computations need only execute simple matrix operations, the computation 
of spin-axis attitude is beset with the burden of solving complex 
relations from spherical trigonometry. 

Since spin-axis attitude is usually not computed frequently, the 
need for efficient algorithms is not immediate, at least not for ground 
support systems. The determination of the spin -axis attitude from batch 
measurements of arc lengths and dihedral angles has become highly 
standardized and reliable^ and there is no obvious need to replace this 
software in normal ground support operations. 

The need for more efficient algorithms lies in two areas: 1) the 

eventual implementation of spin -axis attitude computation in onboard 
microprocessor-based attitude determination systems; and 2) the computa- 
tion of spin -axis attitude accuracies, which imposes a far greater 
computational burden than computing just the attitude due to the greater 
number of terms and because the computation of the attitude covariance 
involves implicitly the computation of derivatives of the attitude. 

The large computational burden imposed by the need to solve 
spherical trigonometric equations in the computation of spin-axis 
attitude covariances is evident in the work of Wertz and Chen,^**^"® 
the most complete and careful work to date. The difficulties which are 
encountered in this approach are of two kinds: 1) the complexity of the 

trigonometric relations, themselves, and 2) the fact that for certain 


16-3 



cases the representation of the quantities being calculated becomes 
indeterminant while the quantities themselves are well defined. This 
last difficulty is simply a manifestation of the fact that the 
representation of rotations by Euler angles is sometimes ambiguous and is 
overcome in the same way, namely, by changing the representation. 

The need for computing spin-axis attitude covariance matrices is 
two-fold. Firstly, it is necessary to be able to assess the accuracy of 
a spin-axis attitude computation during the spacecraft mission. 

Secondly, it is important to be able to predict spin-axis attitude 
accuracies for mission planning, particularly in the determination of 
launch windows. For an example of launch window computations using the 
geometrical approach see Chen.^ 

The purpose of the present work is to develop algorithms for 
computing spin-axis attitude and the associated covariance matrix without 
relying as heavily as do current methods on the solution of trigonometric 
equations. A completely vectorial approach is, of course, not possible 
owing to the nature of the measurements themselves. However, in large 
degree many of the trigonometric equations can be abandoned with the 
result that the spin-axis attitude and, particularly, the covariance 
matrix can be computed more efficiently. 

The types of measurements studied here are of two kinds: 

measurements of arc length , which will always be the angle 
between the observed direction and the spin axis. 

measurements of dihedral angles , i.e., the angle between two 
planes, where the line of intersection is assumed to be the 
spin axis.® 

Dihedral angles, in general, are measured by observing two crossing 
times in the spacecraft and multiplying by the angular velocity. Arc 
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lengths may be measured in a variety of ways, for example, by direct 
sighting (as Of the Sun or a star) or by measuring the component of a 
vector along the spin axis (e.g., the magnetic field vector). The 
measurement of the nadir angle is hybrid in that an arc length (the nadir 
angle) is determined from the measurement of a dihedral angle (the Earth 
width). It is the measurement of the nadir angle which is the source of 
most of the computational complexity. 

Estimation algorithms may be classified either as deterministic 
(usually single-frame, i.e., single-time) algorithms, in which a minimal 
subset of the available data is chosen to compute the spin-axis attitude, 
or as optimal (batch) algorithms, in which a larger quantity of data is 
used from which one computes a "best" result. Three cases are treated in 
this report 

1) A deterministic estimator using two arc -length measurements, 

2) A deterministic estimator using the measurements of two arc 
lengths and the included dihedral angle. (Since in this 
case the spin-axis attitude is over-determined the question 
of optimality is also discussed.) 

3) An optimal batch estimator utilizing any number of 
measurements of dihedral angles and arc lengths. 

In each case the covariance analysis is presented in detail. 

In the appendix the measurement of the nadir angle is presented. It 
is at this point that trigonometric relations cannot be avoided, at least 
in so far as measuring instruments (horizon scanners) are presently 
constructed. The treatment is similar to that of Wertz and his 
collaborators (Ref. 2) but a method is given for avoiding sign 
ambiguities. 
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The treatment of single-axis attitude estimation presented here 
complements that of Wertz. The advantage of Wertz's treatment is that 
the variances along two great circles of the celestrial sphere 
intersecting at the direction of the spin axis and the dihedral angle 
between these two circles (the correlation angle) is given fairly 
directly. Much less direct is determining the covariance of the 
spin-axis vector in inertial space. This part of the calculation falls 
out simply in the present formalism. 

The results presented here are quite simple although they do not 
seem to be generally known. An important result, which is demonstrated 
here, is that little accuracy is lost by relaxing the constraint in the 
optimization that the spin -axis vector be a unit vector and then 
unitizing post hoc. This is responsible for a great deal of 
simplification of the methods presented here, especially for batch 
estimation. 


II. Single-Frame Spin-Axis Estimation from 
the Measurement of Two Arc Lengths 

Consider the simplest case in which the measured quantities are e, 
the Sun angle (the angle between the spin axis and the Sun vector), and 
n, the nadir angle (the direction between the spin axis and the nadir 
vector). The case where one of these measurements is replaced by the 
magnetic field angle is analogous. 

Let ^denote the Sun unit vector, E the nadir vector, and ^ the spin 
axis. Then 

S*n = cos 0 = c- (la) 

E*n = cos n = c_ (lb) 

W bM# I- ' ' 
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The direction of the spin-axis can then be determined simply by using a 
method that has been published recently by Grubin,^ though it has been in 
use since the beginning of the space program and probably has been known 
for several hundred years. 

If ^and J^are not parallel, then it is always possible to write 



i = aJ + a^E + a„ i X E 

(2) 

The problem is 

now to determine the coefficients a^, a^, a^^. 


From Eqs, 

(1) and the normalization condition we have 



C3 = H = + a^d-D 

( 3 a) 


= id= 85 ^E 

( 3 b) 


1 = i.-a = 3s + + 2a3a3(S*E) + a^|i x J|2 

( 3 c) 


which have the solution 




li^i 


2 


■ I *. ^ ?.|2 


S X E 




S X E 


- (Cs - 2C5C^(H) + 


(4a) 


(4b) 


(4c) 
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Note that there are two possible solutions for n. These are shown 
geometrically in Figure 1. 

It will be convenient to define the following quantities 





( 6 ) 


where the tilde below the letter denotes a two-dimensional ‘vector or a 
2x2 matrix, 

Eqs, (4) can now be written 

a = U c ( 7a ) 



S X E 


[1 


- c'^U 


(7b) 


The covariance analysis is now straightforward. Define the three 
vector 



( 8 ) 



Then the covariance matrix of the measurements is given by 

s <6c 6c^> (9) 

where the bracket denotes the expectation value and 6c is the error in c. 
The covariance matrix of the spin-axis direction in the non-orthogonal 
coordinate system is 

Pg = <«a (10) 

and in an orthogonal coordinate system 

( 11 ) 

Substitution of Eqs. (7) in Eq. (10) gives readily 



wi th 

M = <6a 6a^> = U P^ 

~ -c ~ 

V = M b 
S = b^ M b 

16-9 


(13a) 

(13b) 

(13c) 
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where 


B = 


I 


I I b 

I ~ 



(18) 


Equations (17) and (14) may now be combined to give 


where 



2 

I 


i=l j=l 




Xi = 

w-1 


I.+ 


(19) 


(20a) 



» E + 


I) 


Eq. (16) 1s again satisfied since 


(20b) 


,X.-n = 0 i=l,2 (21) 

III. Single-Frame Spin-Axis Estimation from the Measurement 
of Two Arc Lengths and the Included Dihedral Angle 

The ambiguity in determining the spin-axis observed in the previous 
section is removed if the included dihedral angle is also measured. The 
dihedral angle t is defined as the angle between the (S,n) and (£,n) 
planes and is easily shown to be given by 


sin i|> 




■'(i-(5-i)^)(i-(£-a)^; 


(22a) 
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(22b) 


cos = 


( 1 * 1 ) - 

Al-{^n)'^)(l-(E•n)^) 


tan ij> = 


A A A 

Ji* (1 £) 


(H) - {H){|*i) 


(22c) 


The geometry is depicted in Figure 2. 


To determine the spin axis attitude it will be convenient to define 


=-i{l-Cs){l-c|) sin ♦ 


(23) 


and 


c = 


(24) 


The vector a is now determined by four equations 




(25a) 




(25b) 


A 

= s 


X eI^ a 


(25c) 


1 = * a| * 2 ajajCS-l) + ag 


(25d) 
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The three components of ^ are now overdetermined. The most convenient 
solution is obtained by solving the first three equations, which are 
linear, leading to 

a = U c (26) 

WWW 

where 


U = 


S X E 


/N A 

1 0 

-(H) 1 0 

0 0 1 


(27) 


The spin-axis in given by this however, is not properly normalized 
since the measurements are not exact, A properly normalized spin-axis 
vector is then obtained by simply normalizing the solution 


"=ii/|a| 

The covariance matrix of a is given simply by 


(28) 


P = U 
a c 


(29) 


and the covariance matrix for the unnormal ized spin -axis is given by 


P„.TP^tT (30) 


similarly to Eq. (14). The covariance matrix of the properly normalized 
spin-axis vector is recovered simply as 



QPnQ 


(31) 
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where 


Q = I - n (32) 

It is well to ask how good is the approximation of ignoring the 
normalization condition and then normalizing the solution post hoc. 
Instead of this seemingly brutal approach one can find the best solution 
to Eqs. (25abc) subject to the constraint of Eq. (25d), i.e., one seeks 
to minimize the loss function 


L(a) = (c-Aa)"^ (c-Aa) 


subject to the constaint 


where 


a^ A a = 


1 



1 (S*E) 0 

(S-i) 1 0 


(33) 


(34) 


(35) 


The solution is straightforward and yields 

a„pt - (A c (36) 

where X is the Lagrange multiplier for the constraint and from Eq. (34) 
is the root of the equation 

c c 

which yields the smallest value of the loss function. 
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Equation (36) may be rewritten 


Jopt ° (I-lPcU)-'a 


(38) 


where a is given by Eq. (26). Since is expected to be close to a, 
it follows that XP^U must be small. An approximate solution for a^p^ can 
be obtained by expanding Eqs. (37) and (38) in xp and solving. This 
yields 


Now 


-opt 


1 p n a 

2 aVA-^a “ 

c 


(39) 


<l-a^Aa> = Tr(Pj.U) 


(40a) 


<(l-a^Aa)^> = 4 aVa 

' MM' «nf ^ mM 


(40b) 


so that the additional root mean square (rms) error in ^ when optimality 
is not taken into account is of the same order of magnitude as the rms 
error in the cosine measurements. However, the source of this additional 
error, as shown by Eqs. (40) is the error in the normalization. Hence 
this error will be greatly reduced when the unit vector is normalized. 

IV. Batch Estimation 

The value of avoiding trigonometric expressions becomes more obvious 
in dealing with batch estimation. The computational advantage of the 
present approach over the geometrical approach^ is substantial. 
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For batch estimation the non -orthogonal basis cannot be used since 
only the Sun vector is constant (and then only for relatively short data 
spans). The present treatment focuses on the case where the measurements 
consist of two arc lengths and the included dihedral angle. The 
extension to other cases is straightforward. 


Let Cj(i), Cg(i), Cj^(i) be a series of measurements of the Sun 
projection, the nadir projection, and the Sun-nadir dihedral angle, 
respectively. Then the best solution for the spin-axis is obtained by 
minimizing 



subject to the constraint 



( 42 ) 


In order to decrease the number of subscripts in the expressions it has 
been assumed that each data type is available at each time and that each 
measurement type has a single characteristic error. Except for a 
proliferation of subscripts the expressions which follow are not changed 
when this assumption is removed. 
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The minimization of Eq. (41) subject to the constraint is 
straightforward and leads to 

n,= (M-Xir’v (43) 


where 


M = 


'"1 "s "e “n 


(44a) 


V = 


1=1 <>s ‘^E 


(44b) 


and X is the root of 


^1= 1 

(M-XI)2 


which leads to the smallest value of Eq. (41), 


(45) 


As in the previous section it can be expected that the constraint 
can be ignored (X»«0) and the solution be approximated by 


i = i!/|l»| (46) 

where 

n = 1 (47) 
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This approximation has been tested for one spacecraft^® and been observed 
to be quite good. The covariance of j]^ is given by 

P„ = «■' (48) 

and the covariance of the normalized solution is given again by 

P = Q P„ Q (49) 

n 


V. Measurement Errors 

The computation of the spin-axis covariance matrix requires as 
input a model for the covariance matrix of the cosine measurements. 
Expressions are derived here for computing these for the case of Sun and 
Nadir measurements. The treatment when one of these measured quantities 
is the magnetic field is treated in the same way. 

Sun Measurements 

The quantity measured is usually the Sun angle, 3. Hence, 

5Cj = -sinB6& (50) 


Nadir Measurements 


If the spacecraft has angular velocity w, then the Earth width is 
given by 


« = “(tg-tj) (51) 

where tj and t^ are the in- and out-triggering times, respectively, of 
the Earth scan (for a momentum-wheel mounted scanner, o) will be the 
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angular velocity of the momentum wheel). 


Then, using the results from the appendix 


= 6cosn 


3cosn 

H «cos^ 


9C0S^ 


Sinn ^ 

cotY - cotn ^ ? 


= ^ (sin (6t^ - 6t.) (52) 

2 cotY - cotn ' 2^ ' 0 I' ^ ' 

where y is the scan-cone half angle. 

Dihedral Angle Measurements 

The dihedral angle t is determined from the time interval from the 
Sun crossing to the mid-point of the horizon scan 

= 0)[tj - -^(tp + tj)3 (53) 

Thus, (3,Ji,'l») or (0,n,t) is a set of statistically independent 
variables. The "dihedral cosine” Cn^, however, is given by 

C|^| = sin3 Sinn sirnj^ (54) 


hence 


6C|y| = C|nj[cote 60 + cotn 6n + cott 6i|>] (55) 

From Eqs. (50-55) the covariance matrix can easily be calculated. 
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To a large degree, much of the trigonometric complexity which has 
been removed from the attitude solution has simply been shifted to the 
computation of a derived measurement covariance matrix. There is, 
however, a substantial gain because the covariance matrix need not be 
computed to the same degree of accuracy as the spin-axis attitude 
itself. Hence, a great deal of computational approximation is possible, 
such as approximation of the trigonometric functions by simple rational 
functions. 


Appendix - Measurement of the Nadir Angle 

Because the Earth is an extended body the nadir vector is not 
measured directly but determined from measurements of the Earth width. 
Earth widths are measured by a horizon scanner, which effectively is a 
sensor mounted on a rotating cone (of half-cone angle y) about the 
spacecraft spin axis, which detects the crossings of the Earth horizon on 
the scan cone. The Earth has an effective angular radius of p, which is 
a function of altitude and (for a non-spherical Earth) latitude. The 
Earth width is the dihedral angle between the in- and out-crossings (Hj 
and Hq) the horizon by the scanner and is denoted by n. These quantities 
are related by the spherical law of cosines^ 

cosp = COSY cosn + sinY sinn cos(J2/2) (A-1) 

The geometry is depicted in Figure 3. 

Eq. (A-1) may be solved to give 


cosn 


COSP COSY ± 


sinp cos 


(n/2) 4 


-cos P 


A 


(A-2) 
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where 


A = cos^p + sin^Y cos^(n/2) (A-3) 

The sign ambiguity may be eliminated if another measurement is 
present, say that of the Sun angle, 3, and the Sun-Earth dihedral angle, 
Y. Let C be the arc length from the Sun direction to the mid scan point 

cos5 = cos3 COSY + sin3 sinY cosi{> (A-4) 

Then it is possible to show that the underdetermined sign in Eq. (A-2) 
must be the same as that of 

(cos3 - cosy) (E*S - cos5) 

Alternatively, one may consider simultaneously Sun and horizon 
measurements. This leads to three simultaneous equations 

cosB cosn + sin$ sinn cost = (A-5a) 

COSY cosn + sinY sinn cos(n/2) = cosp (A-5b) 

2 2 

COS n + sin n = 1 (A-5c) 

Equation (A-2) was obtained by solving Eqs. (A-5b) and (A-5c) 
simultaneously. One could just as easily solve Eqs. (A-5a) and (A-5b) 
for cosn and sinn. The result will not necessarily satisfy Eq. (A-5c) 
but the two equations have the advantage of being linear. The solutions 
can then be renormalized to satisfy Eq. (A-5c). 

This approach of ignoring the proper normalization for the 
trigonometric functions has another advantage in that a simultaneous 
solution to Eqs. (A-5b) and (A-5c) may not exist in certain extreme cases 
because the measurements are not exact. By solving Eqs. (A-5a) and 
(A-5b) a solution will always exist. 
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There is, however, one clear disadvantage. If Eq. (A-2) is used then 
e, n, and will be statistically independent. If, however, the linear 
equations are solved, n will be correlated with 0 and Thus, the 
simplicity gained in computing cosn is counterbalanced by greater 
complexity in computing the measurement covariance matrix P . 



Figure 1 

Single-Axis Attitude from Two 
Arc-Length Measurements 
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Figure 2 

Single-Axis Attitude from Measurements 
of Two Arc Lengths and One Dihedral Angle 


A 

A 



A 

A 



Figure 3 

Geometry for Nadir-Angle Determination 
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